AD-A042  549  AlR  FORCE  FLIGHT  DYNAMICS  LAB  WRIGhT-PATTERSON  AFB  OHIO  F/6  4/2 

A DESCRIPTION  OF  THE  ATMOSPHERIC  TURBULENCE  ENVIRONMENT  DERIVED— ETC (UJ 
APR  77  P L hasty 

UNCLASSIFIED  AFFDL-TR-77-4  ML 


1 OF  1 
A0AO42540 

a 

It’ 

r 

1 

1 



END 

niMD 

8*77 

DDC  FILE  copy/  ADA042549 


A DESCRIPTION  OF  THE  ATMOSPHERIC  TURBULENCE 
ENVIRONMENT  DERIVED  FROM  THE  CRITICAL 
ATMOSPHERIC  TURBULENCE  (ALLCAT)  PROGRAM 


STRUCTURAL  INTEGRITY  BRANCH 
STRUCTURAL  MECHANICS  DIVISION 


APRIL  1977 


TECHNICAL  REPORT  AFFDL-TR-77-4 

FINAL  REPORT  FOR  PERIOD  APRIL  1963  - DECEMBER'S 


Approved  for  public  release;  distribution  unlimited 


o ^ 

V-'  _ 


5 IQIT 


AIR  FORCE  FLIGHT  DYNAMICS  LABORATORY 
AIR  FORCE  WRIGHT  AERONAUTICAL  LABORATORIES 
AIR  FORCE  SYSTEMS  COMMAND 
WRIGHT-PATTERSON  AIR  FORCE  BASE,  OHIO  45433 


> 


NOTICE 


When  Government  drawings,  specifications,  or  other  data  ^tre  used  for  any  purpose 
other  than  in  connection  with  a definitely  related  Government  procurement  operation, 
the  United  States  Government  thereby  incurs  no  responsibility  nor  any  obligation 
whatsoever;  and  the  fact  that  the  government  may  have  formulated,  furnished,  or  in 
any  way  supplied  the  said  drawings,  specifications , or  other  data,  is  not  to  be 
regarded  by  implication  or  otherwise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights  or  permission  to  manufacture, 
use,  or  sell  any  patented  invention  that  may  in  any  way  be  related  thereto. 


This  report  has  been  reviewed  by  the  Information  Office  (10)  and  is  releasable  to 
the  National  Technical  Information  Service  (NTIS).  At  NTIS,  it  will  be  available  to 
the  general  public,  including  foreign  nations. 


This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


PAUL  L.  HASTY 
Project  Engineer 


FOR  THE  COMMANDER 


HOWARD  L.  FARMER,  Col,  USAF 
Chief,  Structural  -Mechanics  Division 


ROBERT  M.  BADER,  Chief  • 
Structural  Integrity  Br. 
Structural  Mechanics  Div. 


Copies  of  this  report  should  not  be  returned  unless  return  is  required  by  security 
considerations,  contractual  obligations,  or  notice  on  a specific  document. 

AiR  FORCE  - 28  JUL  77  - 500 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Date  Entered) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


I.  REPORT  NUMBER 


! iiJ- 


AFFDL-TR-77-4 


2,  GOVT  ACCESSION  NO 


(l<^ 


E.  T IT  <«***'<*>' 


lA  DESCRIPTION  OF  THE  ATMOSPHERIC  TURBULENCE  ENVI-  , 
RONMENT  DERIVED  FROM  THE  CRITICAL  ATMOSPHERIC  / 
TURBULENCE  (ALLCAT)  PROGRAM- 


pient's  catalog  number 


flE  MBjjgT  k PERIOPTOllEREi}  _ 

Final  yFepart,  — ‘‘K 

Aprs  )<963  — DecMfir  il672^ 


b.  PE'N'IgRMINU  URB.'HEMUTmW'OWBER'^ 


7.  AUTHOR^*; 

Paul  L.^Hasty 


CONTRACT  OR  grant  NUMBERfaJ 

N/A 


9.  PERFORMING  organization  NAME  AND  ADDRESS 

Structural  Integrity  Branch  (FBE) 

Air  Force  Flight  Dynamics  Laboratory 
Wright-Patterson  AFB,  Ohio  45433 


10,  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  & WORK -«^|T  NUMBERS^.  _ _ _ 

62201F/TH^01/13 


I I.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Air  Force  Flight  Dynamics  Laboratory  I' 

Wright-Patterson  Air  Force  Base,  Ohio  45433 


14.  MONITORING  AGENCY  NAME  ft  AODRESSC//  diUerent  {rom  Controlling  Otftce) 


1 2 -RgPOT^T’OAT'C  ' 

AprB  B77- 

13.  NUMBER  OF 

'fifl  ^ 


15.  SECURITY  CLASS 

UNCLASSIFIELD 


I5a.  DECLASSIFICATION  'DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (of  this  Report) 

Approved  for  public  release;  distribution  unlimited 


o 


17.  DISTRIBUTION  STATEMENT  (of  fhe  abstract  enferad  In  Block  20,  If  d//f«rent  from 
\ 

G 

19.  KEY  WORDS  (Continue  on  reverse  side  if  necessary  and  identify  by  block  number) 

Turbulence  Gust  Velocities 

Gusts 

Atmospheric  Environment 
Power  Spectral  Densities 


2e?r  ABSTRACT  (Continue  on  reverse  side  If  necessary  and  Identify  by  block  number) 

This  report  summarizes  the  important  findings  of  the  individual  projects  of 
the  Critical  Atmospheric  Turbulence  (ALLCAT)  Program.  The  basic  gust  parameters 
are  presented  for  the  various  altitude  regions.  Power  spectral  densities  are 
presented  and  the  results  of  the  investigation  are  presented  in  a form  useful 
for  gust  design  procedures.  Scale  lengths,  turbulence  intensity,  and  percentage 
of  flight  time  in  turbulence  are  presented  as  a function  of  altitude. 


DD  1473  EDITION  OF  I NOV  6S  IS  OBSOLETE 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  flFh-n  Pmlm  FnIrrrdI 


f 


AFFDL-TR-77-4 


FOREWORD 

This  report  was  prepared  by  Mr.  Paul  L.  Hasty,  Aerospace  Engineer 
in  the  Loads  and  Response  Prediction  Group,  Structural  Integrity 
Branch  (FBE),  of  the  Structural  Mechanics  Division,  Air  Force  Flight 
Dynamics  Laboratory  at  Wright-Patterson  Air  Force  Base,  Ohio.  The  work 
was  performed  under  Project  1367  "Structural  Integrity  for  Military 
Aerospace  Vehicles",  Task  136701  "Structural  Flight  Loads  Data". 

The  Critical  Atmospheric  Turbulence  (ALLCAT)  Program  was  authorized 
and  funded  mostly  under  the  Advanced  Development  Program  683E.  The 
research  was  conducted  from  April  1963  through  December  1972. 

This  report  presents  a description  of  gust  data  obtained  from 
flight  test  programs.  The  results  are  presented  in  a form  useful  for  the 
design  of  aerospace  vehicles.  This  report  is  a concise  summary  of  the 
individual  projects  conducted  under  the  ALLCAT  Program  or  programs 
directly  related  to  it. 
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SECTION  I 
INTRODUCTION 

The  Air  Force  Flight  Dynamics  Laboratory  conducted  a large  scale 
program,  the  Critical  Air  Turbulence  (ALLCAT)  Program,  to  develop  a 
valid  model  for  the  description  of  the  turbulence  environment.  Large 
quantities  of  data  have  been  obtained  for  developing  improved  gust 
criteria,  and  updating  turbulence  design  criteria.  This  report  presents 
the  significant  results  of  this  turbulence  investigation  and  related 
geophysical,  meterological , and  climatic  conditions.  Values  of  suggested 
basic  turbulence  parameters,  required  for  design,  are  presented.  The 
plan  developed  to  define  models  of  atmospheric  turbulence  over  a wide 
range  of  altitudes  was  as  follows: 

I 
i 

Program 

Take-Off  and  Landing  Critical  Atmospheric 
Turbulence  (TOLCAT) 

Low  Level  Critical  Atmospheric 
Turbulence  (LO-LOCAT) 

Medium  Altitude  Critical  Atmospheric 
Turbulence  (MEDCAT) 

High  Altitude  Critical  Atmospheric 
Turbulence  (HICAT) 


4 


4 


Altitude  (ft)  Above  Ground 
0 - 250 


3 


250  - 1 000 

20.000  - 40,000 

40.000  - 70,000 
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SECTION  II 
PROJECT  TOLCAT 

Turbulence  encountered  during  take-off  and  landing  can  present 
serious  hazards  to  aircraft.  Both  temporal  and  spatial  descriptions  of 
turbulence  are  required  for  design  problems.  The  Take-Off  and  Landing 
Critical  Atmosphereic  Turbulence  Project  (TOLCAT)  was  initiated  to 
develop  experimentally  the  necessary  knowledge  of  turbulence  in  the 
atmospheric  boundary  layer,  its  dependence  on  height,  terrain  features, 
meteorological  conditions,  and  to  determine  the  interrelationship  between 
temporal  and  spatial  characteristics. 

An  experimental  investigation  was  conducted  to  develop  optimum 
methods  and  techniques  for  data  acquisition  and  analysis  required  to 
define  the  atmospheric  boundary  layer  turbulence  environment.  The 
investigation  was  specifically  designed  for  V/STOL  aircraft  flight 
control  and  stability  studies  and  to  develop  valid  structural  design 
criteria.  Measurements  of  turbulence  were  made  from  various  arrays  of 
sensors  mounted  on  towers  to  demonstrate  methods  of  describing  the 
temporal  and  spatial  characteristics  of  turbulence.  The  towers  were 
located  near  Hanford,  Washington.  Measurement  and  data  reduction 
techniques  for  sonic  and  three-propeller  anemometers  were  optimized  to 
assure  accurate  measurements  of  the  vertical  and  two  horizontal  wind 
components  from  the  sensors.  Table  1 summarizes  the  experimental 
arrays. 

Experiments  were  conducted  from  a T-shaped  array  consisting  of  four 
sensors  mounted  on  booms  to  simulate  two  wing  positions  (9  meters  apart), 
fuselage  position,  and  tail  position  (6  meters  from  wings)  on  an  aircraft 
and  thereby  compute  correlations  for  various  gust  and  shear  combinations 
to  which  an  aircraft  would  respond. 

On  another  tower,  six  sensors  were  spaced  logarithmically  along 
a vertical  line.  Sonic  anemometers  were  installed  at  60,  15,  and  7.5 
meters,  and  Gil  1 anemometers  were  installed  at  30,  3.75,  and  1.88  meters 
above  the  ground. 
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TABLE  1 

EXPERIMENTAL  ARRAYS 


Number 

Number 

of 

Towers 

Height 

(Meters) 

Spacing 

1 

4 

58 

Horizontal,  four  sensors,  one  on 
each  tower  which  are  equally  spaced 
(223  meters) 

2 

1 

58 

Horizontal,  T-shaped,  two  sensors 
9 meters  apart  and  two  sensors 
6 meters  apart  - simulate  aircraft 

3 

8 

4 

Horizontal,  eight  sensors,  log  spaced 
over  total  of  1023  meters 

4 

1 

12 

Horizontal,  six  sensors,  1.5  meters 
apart 

5 

1 

1.88  to  60 

Vertical,  six  sensors,  log  spaced 

6 

9 

15 

Horizontal,  L-shaped,  seven  sensors 
log  spaced  on  principal  leg  (128 
meters)  and  sensors  at  12  and  60 
meters  on  short  leg 

7 

9 

30 

Same  as  Array  6 

8 

1 

131 

Single  sensor 
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An  L-shaped  array  included  nine  sensors  horizontally  dispersed.  The 
principal  leg  of  the  L consisted  of  seven  logarithmically  spaced  positions 
with  separation  distances  between  consecutive  sensors  varying  from  4 to 
128  meters.  The  purpose  of  this  line  of  sensors  was  to  make  possible 
the  comparison  of  time  and  space  analyses  when  the  wind  blew  from  various 
azimuth  angles  to  the  principal  line.  The  short  leg  of  the  L-array  had 
two  positions,  12  and  60  meters  from  the  origin.  The  purpose  of  the 
short  line  of  sensors  was  to  make  possible  comparison  of  correlations 
in  two  directions  at  12  and  60  meters.  The  same  tower  locations  were 
used  for  both  the  15-meter  and  the  30-meter  L-arrays.  The  only  difference 
was  the  heights  at  which  the  sensors  were  mounted.  The  purpose  was  again 
to  obtain  spatial  correlations  and  cross  spectra  to  compare  with  time 
series  analysis. 

In  the  final  experiment,  a single  sensor  was  mounted  atop  the 
125-meter  tower  in  order  to  obtain  samples  of  data  at  elevations  above 
the  surface  layer,  where  the  existing  turbulence  is  related  less  to 
the  local  wind  shear  and  more  to  convection  and  vertical  transport  of 
turbulence  from  below. 

Several  tests  were  conducted  with  the  sensors  in  each  of  the  arrays 
described  above.  A list  of  twelve  data  segments  chosen  for  analysis 
on  the  basis  of  their  high  quality  is  given  in  Table  2.  In  the  first 
column,  the  first  digit  gives  the  array  number  (from  Table  1)  and 
second  digit  gives  a serial  number  identifying  the  particular  test 
within  the  series. 

Both  analog  and  field  digital  magnetic  tape  recording  options  were 
demonstrated  to  provide  flexibility  in  measurement  array  conf igurations 
and  to  optimize  recording  capabilities  for  a variety  of  measurement 
requirements.  Data  reduction,  processing,  and  analysis  techniques 
were  developed.  The  detailed  results  are  presented  in  References  1 
and  2. 
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TABLE  2 

TESTS  ANALYZED  FOR  TOLCAT 


Mean  Wind 


Test  No. 

Date 

Time  Segment 

Stabi 1 1 ty 

Direction  (Deg) 

Speed  (mps) 

Fleighl  (m' 

15 

681202 

1345-1426 

SI ightly  Stable 

212 

8 

58 

22 

690925 

1552-1647 

SI ightly  Stable 

296 

9 

58 

51a 

700720 

1537-1605 

Unstable 

320 

9 

60* 

51h 

700720 

1619-1647 

Unstable 

304 

10 

60* 

56 

700723 

2143-2238 

Stabi  e 

301 

12 

60* 

58 

700728 

1910-2005 

Neutral 

303 

9 

60* 

63 

710119 

1323-1  351 

SI  ightly  Stable 

300 

6 

15 

71 

710303 

1316-1344 

Unstable 

243 

14 

30 

73 

710324 

0945-1013 

Unstable 

303 

7 

30 

74 

710326 

1013-1059 

Unstable 

200 

13 

30 

83 

710607 

1143-1224 

Unstable 

297 

5 

131 

85 

710621 

1059-1154 

Unstable 

50 

3.6 

131 

*Turbulence  ineasurements  were  taken  at  lower  heights 
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Some  experimental  data  were  collected  from  tower  arrays  as 
shown  i n Tabl e 1 . 

Only  some  general  conclusions  can  be  stated  as  a result  of  the 
TOLCAT  investigation.  Initial  studies  of  Taylor's  hypothesis  showed 
that  autocorrelation  functions  scaled  by  the  mean  wind  speed  to  the 
expected  spatial  functions  agreed  quite  well  with  the  measured  spatial 
correlations.  This  hypothesis  states  that  the  space  and  time  auto- 
correlations (or  equivalently  spectra)  are  equal  (with  the  transformation 
x^  = Ut).  U is  the  speed  of  the  airplane  with  respect  to  the  atmosphere, 
and  x^  refers  to  the  three  space  directions.  The  use  of  the  hypothesis 
allows  one  to  obtain  information  by  using  one  probe  and  doing  analysis 
in  time,  instead  of  using  two  probes  and  doing  analysis  in  space.  The 
translation  of  time  series  spectra  to  that  experienced  by  an  aircraft 
holds  only  for  the  case  where  the  aircraft  flies  in  the  direction  of 
the  wind. 

Modeling  of  turbulence,  appropriate  to  take-off  and  landing 
problems,  will  require  an  extensive  measurement  and  analysis  effort  to 
follow  this  development  phase. 
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SECTION  III 
PROJECT  LO-LOCAT 

The  LO-LOCAT  Program  was  established  to  determine  the  turbulence 
environment  below  1000  feet  above  the  ground  utilizing  statistically 
representative  samples  of  turbulence  data  obtained  over  a wide  range  of 
meteor ological , topographical , seasonal,  and  time-of-day  conditions. 

The  following  major  areas  were  studied:  probabilities  of  encountering 
given  levels  of  turbulence  intensity,  frequency  analysis  of  atmospheric 
turbulence  utilizing  the  power  spectra  method,  correlation  of  atmospheric 
turbulence  with  meteorological  and  geophysical  conditions,  and  associated 
forecasting  techniques.  The  basic  information  concerning  the  three 
phases  of  the  LO-LOCAT  Program  are  presented  in  Table  3 from  Reference  3. 

The  LO-LOCAT  data  consist  exclusively  of  clear  air  data.  No  storm 
flights  are  included.  During  the  program,  specific  routes  were  flown 
on  a scheduled  basis  without  regard  to  the  amount  of  turbulence  predicted. 
The  data  sample  can  be  considered  a random  sample  of  the  gust  environment 
of  the  atmosphere  during  nonstorm  conditions.  The  model  used  to 
represent  atmospheric  turbulence  is  considered  to  be  made  up  of  discrete 
patches  of  turbulence  of  different  intensity,  where  each  patch  is  assumed 
to  be  Gaussian  and  stationary.  If  a single  stationary,  Gaussian  patch 
is  considered,  the  average  number  of  gust  peaks  per  unit  distance  which 
exceed  a given  level  of  vertical  velocity  is  given  by 


N(w)  = N e 
0 


-('S) 

2o 


(This  is  true  for  one  particular  turbulence  spectrum).  is  the 
average  number  of  times  per  unit  distance  that  the  gust  velocity 
function  crosses  the  value  zero  with  positive  slope. 
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TABLE  3 

BASIC  LO-LOCAT  INFORMATION 


b 


Phases  I & 

II 

Phase  III 

Aircraft  used  to  obtain  data 

4 C-131's 

1 T-33 

Flight  period 

9-15-66  - 

12-20-67 

8-16-68  - 6-30-69 

1 Number  of  turbulence  data 

1 flights 

1 ,244 

299 

1 Low  altitude  turbulence  data 

1 samples 

8,871 

1,938 

j BREN  tower  samples 

0 

8 

I Thunderstorm  turbulence 

samples 

0 

55 

Wake  turbulence  samples 

0 

46 

Turbulence  sample  average 
length 

21  mi . 

32  mi . 

Turbulence  sample  average 

time  interval 

5.5  min. 

4.5  min. 

Average  turbulence  PSD  wave- 
length band 

33  - 7,980 

feet 

63  - 15,100  feet 

Aircraft  average  ground 
speed 

332  fps 

630  fps 

u time  series  samples  anal- 
yzed statistical ly 

6,553 

1 ,730 

V time  series  samples  anal- 
' yzed  statistically 

6,270 

1,709 

w time  series  samples  anal- 
yzed statistically 

6,508 

< 

1,716 

' u non-contour  samples 

1,113 

0 

V non-contour  samples 

1 ,091 

0 

w non-contour  samples 

1,118 

0 

u rms  values  analyzed 
statistical  ly 

7,670 

1,762 
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TABLE  3 (Contd) 


Phases  I & II 

Phase  III 

V rms  values  analyzed 
statistical ly 

6,621 

1 ,740 

w rms  values  analyzed 
statistically 

7,630 

1,746 

Gust  velocity  samples  for 
which  spectra  were  calculated 

1,775 

904 

Homogeneous  u spectra  analyzed 

1 ,272 

591 

Homogeneous  v spectra  analyzed 

1,061 

524 

Homogeneous  w spectra  analyzed 

1,259 

496 

u spectra  obtained  during  non- 
contour flight 

207 

0 

V spectra  obtained  during  non- 
contour fl ight 

200 

0 

w spectra  obtained  during  non- 
contour flight 

210 

0 

Gust  velocity  time  series 
samples  per  turbulence  sample 

33,000 

27,000 

Time  interval  - spectra 
calculations 

0.01  sec. 

0.02  sec. 

Gust  velocity  samples  in 
spectra  calculations 

33,000 

13,500 

Frequency  interval  - spectra 
calculations 

0.0416  cps 

0.0416  cps 

Nyquist  frequency 

50  cps 

25  cps 

Number  of  frequency  estimates 
PSD 

1,200 

600 

Miles  covered  recording 
acceptable  turbulence  data 

135,900 

55,750 

Hours  of  acceptable  turbu- 
lence data 

600 

130 

Spectrum  degrees  of  freedom  - 
individual  samples 

55 

45 

i 

i 


1 
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For  the  complete  model,  the  equations  for  N(w)  become 


N(w) 


2o?  2or 


d L ^1 


= N, 


w 


d^  N^e 


” ? 2 

2o^  20^ 

I^P-je  + P20  + 


d is  the  total  distance  traveled  by  an  aircraft  in  turbulence  and  d^  is 
the  distance  traveled  in  turbulence  of  intensity  o-j  . For  the  limiting 
case  of  a continuous  variation  in  the  rms  value  o,  this  equation  becomes 

, ' 2o2 

N(w)  = / P(o)e  do 

0 

where  p(o)  = probability  density  of  o 

In  a study  performed  by  the  University  of  Dayton  Research  Institute 
(Reference  4),  it  was  found  that  the  cumulative  distribution  of  peak 
count  exceedance  F(x)  is  reasonably  well  represented  by 


-x/b,  -x/bp 

F(x)  = P^e  ' + p2e 


and  the  density  function  by 


f(x) 


'X/b^ 


-x/b2 


For  the  LO-LOCAT  data,  P^  and  P2  represent  the  percentage  of  total 
flight  time  spent  in  nonsevere  and  severe  turbulence,  b^  and  b2 
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are  turbulence  scale  parameters  for  the  individual  probability  dis- 
tributions of  a^.  Using  the  equation 


F(x) 


-x/b,  -x/b- 

P^e  ‘ + P^e 


the  LnF(x)  was  represented  by  a quadratic  of  the  form 
LnF(x)  = A + Bx  + Cx^ 


This  transformation  requires  a representation  of  P^,  P^,  b^  and  b2  in 
terms  of  A,  B,  and  C.  If  C is  negative  this  indicates  that  the 
cumulative  probability  function  is  concave  down.  The  best  fit  that  can 
be  made  to  such  data  is  to  define  P^  = 1 and  P^  = 0.  The  LnF(x)  reduces 
to  the  form:  LnF(x)  = A-mx,  where  m = l/b.|.  When  C is  positive, 

P.| , P^,  b^  and  b2  are  determined  from  the  coefficients  of  the  equation: 

LnF(x)  = A + Bx  + Cx^ 

b.|  is  defined  as  the  negative  reciprocal  of  the  slope  of  LnF(x)  = 

B + 2 CX;  at  X = 0,  the  slope  = B and  b^=  1/B.  To  determine  b2>  the 
negative  reciprocal  of  the  slope  of  LnF(x),  the  point  of  transition 
between  the  influences  of  the  first  and  second  terms  of  Equation 

-x/b,  -x/b« 

F(x)  = P^e  ' + P20 

is  chosen  as  x^  = x^^^/2  where  x^^  is  the  maximum  gust  velocity  encountered. 

1 1 

B + 3 Cx  ■ ‘B  + 3 Cx 
c m 

2 

B . 3^ 
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Recommendations,  from  Reference  4,  are  given  for  Pp  P^,  and  the 

von  Karman  L for  the  0-1000  ft  altitude  region  in  nonstorm  flight 
conditions 


PlSxp  - + P^exp  - g-  = exp(Bx^  + Cx^ 


The  section  gives  recommendations  for  P^ , P^,  b-| , b2  and  the 
von  Karman  L for  the  0-1000  ft  altitude  region  in  nonstorm  flight 
conditions  (Reference  4).  Also,  probability  of  exceedance  curves  are 
presented  for  various  environmental  conditions.  Figure  1 shows  the 
cumulative  exceedance  curves  for  the  Phases  I and  II  and  Phase  III  data. 
For  Phases  I and  II  a maximum  gust  of  66  ft/sec  in  the  vertical  component 
was  observed  during  a noncontour  flight.  During  Phase  II  the  maximum 
gust  was  76  ft/sec  over  high  mountains.  The  cumulative  exceedance  curves 
for  all  Phases  I and  II  flight  legs  at  250  ft  and  for  all  legs  at  750  ft 
are  presented  in  Figure  2.  Little  difference  is  seen  between  the  two 
curves.  The  cumulative  probability  curves  of  each  component  of  Phase  1 
and  Phase  II  data  is  shown  in  Figure  3.  Similar  results  were  obtained  for 
the  Phase  III  data 


The  cumulative  probability  exceedance  curves  for  all  Phase  III 
flight  legs  have  been  shown  to  be  considerably  larger  than  those  of 
Phases  I and  II.  If  one  considers  only  those  Phase  III  legs  not  flown 
over  high  mountains,  there  is  good  agreement  with  Phases  I and  II. 

Hence,  it  was  decided  to  combine  Phase  III  data,  excluding  high  mountain 
data,  with  all  of  Phases  I and  II.  Figure  4 shows  the  cumulative 
probability  curve  for  all  of  Phases  I,  II  and  III,  excluding  high  mountain 
data.  The  corresponding  P's  and  b's  are  listed;  Phase  III  high  mountain 
data  are  also  shown.  A decomposition  of  the  data  of  Figure  4 into  legs 
flown  in  very  stable  air  and  those  not  flown  in  very  stable  air  is  given 
in  Figure  5.  Table  4 summarizes  the  recommended  P’s  and  b's  for  the 
0-1000  ft  altitude  range.  The  exceedance  curves  of  Figure  4 have  been 
presented  by  longitudinal,  lateral,  and  vertical  components  in 
Figures  6 and  7. 
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2 6 10  14  18  22  26  30  34  38  4 2 46  50 

Gust  Velocity  ft/sec 


Figure  3.  Cumulative  Probability  of  Exceedance  Phases  I and  !I 
Longitudinal,  Lateral,  and  Vertical  Components 
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Al  I Phases  I and  1 1 
and  Phase  III  Not 
High  Mountains 


6 10  14  18  2 2 26  30  34  38  4 2 4 6 50  54  58  62  66  70  74  78 

Gust  Velocity  ft/sec 


Figure  4.  Ciiiiiuidtive  Probability  of  Exceedance  for  Phases  I and  II 
and  Not  High  Mountains  Phase  III  Legs;  Phase  III  High 
Mountains  Legs 
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! 


All  Components  Combined 
Phases  I and  II  ond  Phase  III  Not  High 


Gust  Velocity  ft/sec 


Figure  5.  Cuinuliitive  Probability  of  txceedance,  Very  btablp-vs-Nnt 
Very  Stable  Legs 
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TABLE  4 

RECOMMENDED  P's  AND  b's  FOR  0-1000  FOOT  ALTITUDE 
NONSTORM  CONDITION 

Data  Includes  All  Phases  I and  II  Legs 
and  Phase  III  Legs  Not  Over  High  Mountains 

All  Components  Combined 


f’l 

P2 

^1 

^2 

All  Legs 

0.99986 

0.00014 

2.605 

8.684 

Very  Stable 

0.99750 

0.00250 

2.099 

5.211 

Not  Very  Stable 

0.99978 

0.00022 

2.678 

8.033 

Data 

Includes  Phase  III  High  Mountain 
All  Components  Combined 

Legs  Only 

Pi 

P2 

h 

^2 

High  Mountains 

0.85 

0.15 

3.365 

6.947 

High  Mountains 
Very  Stable 

0.983 

0.017 

2.931 

8.033 

High  Mountains 
Not  Very  Stable 

0.88 

0.12 

3.387 

7.165 
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Phase  III  High  Mountame 


Longitudinal 


Pi  = 0.88 
Pp  =0  12 


P,  =0  88 


P,  =0  30 


bj  = 6.  947 


b,  = 3 995 

bj  = 7 186 


b|  = 3 995 
b,  = 7 186 


Longi  tud  ina  I 


2 6 to  14  18  22  26  30  34  38  42  46  50  54  58  62  66  70  74  78 

Gust  Veloci  ty  ft/sec 


Figure  7.  Cumulative  Probability  of  Exceedance  for  Phase  III  High 
Mountains  Legs  by  Component 
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The  von  Karman  L's  were 
the  results  are  presented  in 


obtained  for  all  components  of  the  data  and 
Table  5. 


A complete  record  of  the  LO-LOCAT  Program  is  presented  in 
References  3 through  7. 


TABLE  5 

MEAN  VON  KARMAN  L's 
°hases  I and  II 


Component 

A1 ti tude 
250  ft.  750 

All  Components 

342 

442 

Longitudinal  Component 

368 

411 

Lateral  Component 

361 

465 

Vertical  Component 

302 

449 

Phase  III 


Component 

Terra i n 
High  Mountains  Not 

High  Mountains 

All  Components 

704 

564 

Longitudinal  Component 

813 

687 

Lateral  Component 

704 

555 

Vertical  Component 

579 

430 

1 


J 

1 

■J 

1 

A 


\ 
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SECTION  IV 
PROJECT  MEDCAT 

The  objective  of  the  MEDCAT  Project  was  to  measure  and  characterize 
the  turbulence  in  the  20,000-40,000  ft  regime.  Flights  were  planned  to 
penetrate  those  regions  predicted  as  being  most  likely  to  contain 
turbulence.  Hence  the  percentage  of  time  in  turbulence  experienced 
during  this  program  is  not  a valid  estimate  of  an  expected  time  in 
turbulence  during  operational  missions. 

F-106  and  F-lOO  aircraft  flew  a total  of  284  flights  which  resulted 
in  121  hours  and  288  hours  respectively,  of  flight  time  in  search  for 
turbulence.  A complete  description  of  this  program  and  an  analysis  of 
the  data  collected  is  given  in  Reference  8.  The  primary  characterizations 
of  the  turbulence  encountered  are  expressed  in  terms  of  power  spectral 
densities  and  counts  of  true  gust  velocities  of  the  three  turbulence 
components.  The  wave  lengths  represented  in  the  power  spectral 
densities  have  a lower  limit  of  100  ft  (approximately  8 cps)  and  a 
maximum  limit  of  1/20  of  the  length  of  the  processed  run. 

After  the  quality  evaluation  of  the  data,  none  of  the  F-106  flights 
were  processed.  Only  25  flights  were  processed  from  the  F-lOO  data  bank, 
those  judged  to  contain  turbulence  intensity  yielding  rms  values  greater 
than  1.5  ft/sec.  Any  conclusions  concerning  the  character  of  turbulence 
from  the  MEDCAT  data  sample  must  be  made  with  reservations.  Accuracy 
limitations  precluded  the  possibility  of  removing  all  aircraft  motions 
from  the  time  gust  computations,  particularly  that  characteristic  ratio 
of  the  F-lOO  aircraft  due  to  the  short  period  longitudinal  and  dutcli 
roll  modes.  However,  conclusions  of  a general  nature  can  be  drawn. 

Estimates  of  the  following  basic  parameters  are  required  from  the 

collected  data  to  use  the  theoretical  models  which  describe  power  spectral 

5 densities  of  atmospheric  turbulence:  the  slope,  m (of  a log-log  plot 

I of  the  power  spectral  density)  in  the  short  wave-length  region;  the 

scale  of  turbulence,  L;  and  the  mean  square  value  of  the  true  qust 
2 

velocities,  o . from  an  analysis  of  the  MEDCAT  data  it  was  concluded 
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that  the  scale  of  turbulence  L is  > 2000  ft,  m = -1.83  for  the  vertical 
component,  and  the  most  probable  o values  were  4.5  ft/sec.  for  the 
lateral  component  and  3.5  ft/sec.  for  the  vertical. 

In  the  F-lOO  data  sample,  less  than  1 of  the  time  was  in  turbulence 
for  which  at  least  one  component  had  an  rms  value  greater  than  1.5  ft/sec. 
The  data  also  indicate  that  the  turbulence  is  not  isotropic. 

Typical  power  spectral  densities  (PSD)  for  the  lateral  and  vertical 
components  have  shapes  represented  by  the  curves  of  Figure  8.  The 
aircraft  motion  distortions  are  not  as  pronounced  in  the  vertical 
component  as  in  the  lateral  but  on  some  of  the  vertical  spectra  peaks 
and  valleys  quite  similar  to  that  of  the  lateral  component  were  present. 

The  plots  of  the  power  spectral  densities  of  true  gust  velocity 
components  on  log-log  paper  typically  exhibit  a straight  line  portion 
in  the  short  wave  length  region  of  the  curves.  The  slope  of  this  straight 
line  portion  provides  an  estimate  of  the  exponent  from  assumed  turbulence 
models  in  this  inertial  subrange  of  the  power  spectral  density.  There  was 
considerable  variation  in  slope  between  the  individual  runs.  The  slopes 
of  the  lateral  component  ranged  from  -1.2  to  -2.2  with  a mean  of  -1.77, 
while  those  of  the  vertical  component  ranged  from  -1.0  to  -1.9  with  a 
mean  of  -1.41.  Histograms  of  these  slopes  are  presented  in  Figure  9. 

In  order  to  arrive  at  one  spectrum  for  use  in  model  development 

and  design  criteria,  average  power  spectral  densities  were  computed  for 

the  lateral  and  vertical  components.  Each  power  spectral  density  was 

2 2 

first  normalized  by  dividing  the  abscissa  by  ^2000’  *^2000 

associated  with  cut-off  wavelength  of  2000  feet.  A weighted  average  of 
these  normalized  spectra  was  then  obtained,  the  weighting  factors  being 
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Loteroi  Component 


Figure  9.  Histogram  of  Spectral  Slopes 
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dependent  on  the  run  length.  The  average  spectrum  was  computed  at 
approximately  25  equally  spaced  log  k values  by  means  of  the  formula 

n • 


1 1 • 
J 

E 

i = l 


where 

({)(kj)  = value  of  average  normalized  power  spectral  density 
T^.  = length  in  seconds  of  run  i 

'l>j(kj)  = value  of  normalized  power  spectral  density  of  run  i at  kj 

n.  = number  of  power  spectral  densities  for  which  an  estimate 
^ is  available  at  \j 

Averages  were  computed  for  all  for  which  Oj  > 4.  Figure  10  presents 
the  average  normalized  spectra  for  all  lateral  and  vertical  runs  which 
met  editing  criteria. 

Average  spectra  were  also  obtained  for  the  two  components  for  the 
three  altitude  bands,  h < 30,000,  30,000  < h < 35,000,  and  h > 35,000  feet. 
Figures  11  and  12  present  these  average  spectra  for  lateral  and  vertical 
components  respectively. 

The  peak  count  cumulative  distribution  for  all  runs  meeting  the 
editing  criteria  were  processed  into  composite  peak  count  exceedance 
rates  for  each  component  direction,  u,  v,  and  w.  The  resultant  cumu- 
lative peak  count  rates  are  shown  in  Figure  13.  From  the  cumulative 
rate  distributions  the  following  conclusions  were  made: 

(a)  Clear  air  turbulence  is  composed  of  more  than  one  rms 
intensity  for  each  component. 


I 
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Figure  10.  Average  Power  Spectral  Densities;  Normalized  by 
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Average  Normalized  Power  Spectrol  Densities -(cpf) 


Cumulative  Peaks  Per  Second 
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GusI  Velocity  Amplitude  - tps 


Fiiiure  It.  Componont  Effects  on  Gust  Velocity  Peaks 
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b ’bp  , urrence  rate  of  the  longitudinal  components  appears 
to  he  • ■ 'jtit  percentage  of  the  vertical  component  except  at  the 

eAtrtH  e ^ ■ .r  tne  amplitude  scale. 

;■  In  the  very  low  amplitude  regions,  < 2 feet  per  second,  it 
appears  that  some  noise  exists.  These  data  were  not  altered  to  correct 
for  this  noise. 

The  cumulative  distributions  of  the  observed  peaks  are  compared 
to  the  distribution  of  peaks  as  specified  in  Specification  MIL-A- 
8866  (USAF)  in  Figure  14.  The  average  of  the  two  sets  of  values 
and  b^.  for  the  altitude  bands  which  span  the  MEDCAT  regime  were  used. 
Thus,  the  MlL-A-8866  (USAF)  cumulative  distribution  was  determined  by 
means  of  the  formula; 


i f -x/b,  -x/b„' 

F(x)  = p~  + ^2® 

= j^0.0625e‘^^^-‘^^  t 0.0009e"^^^^ 

Figure  14  indicates  that  the  MIL-A-8866  (USAF)  spectrum  is  significantly 
more  severe  than  that  obtained  from  the  MEDCAT  data  for  all  three 
components . 

The  commonly  used  theoretical  results  of  Rice  predict  the  distri- 
bution of  level  crossings  of  a stationary  Gaussian  random  process. 

However,  level  crossing  data  are  not  of  prime  interest  in  design  criteria. 
Thus,  it  is  of  interest  to  compare  the  results  of  level  crossing 
analysis  with  the  results  obtained  by  peak  counting  and  amplitude 
counting.  In  order  to  make  comparison  of  results  meaningful,  the  results 
were  converted  into  the  format  of  cumulative  probability  versus 
amplitude.  Plots  are  presented  for  each  of  the  gust  velocity  components 
using  the  results  from  the  three  counting  methods  in  Figures  15 
through  17. 
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Fiqure  14.  MF.DCAT  Lxcppditnco  Peak  Prnhability  vs.  MIL-A-R866(USAF) 
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Figure  15.  Lonqi  tudi  iid  1 Gust  Cumulative  Probability  Distributions 
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Figure  16.  l.itterdl  Gust  Cumulative  Probability  Distributions 
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riqure  17.  VprticrfI  Gi;?t  C imiu I d t i vo  Frobdhilit.v'  Distritiutions 
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A difficulty  arises  in  the  evaluation  of  observed  rms  values  to 
that  present  in  the  evaluation  of  the  scale  of  turbulence.  The  rms 
values  desired  from  the  time  histories  were  not  comparable  since  the 
individual  runs  were  high  pass  filtered  at  different  cutoff  frequencies. 
It  was  necessary  to  use  truncated  rms  values.  The  distribution  of  these 
values  did  not  adequately  describe  the  parameter  of  the  turbulence 
models  as  these  were  stated  in  terms  of  a (not  a^] . The  truncated  rms 
values  must  be  transformed  into  total  rms  values  by  accepting  a family 
of  theoretical  models  and  determining  the  ratio  of  total  rms  to  truncated 
rms  as  a function  of  the  model.  The  investigators  used  the  Taylor-Bul len 
theoretical  models  to  represent  the  spectra  rather  than  the  Dryden  or 
von  Karman  models.  For  the  Taylor-Bul len  Transverse  Family,  this  can  be 
accomplished  by: 


and 


o^k. 


1 

/ 


(t)*(k,L,m)  = oP(L,m,k, ) 


(j)(k)  = 2a^L 


ri  + 2 (n  + 1)  b^  (2tt  kp^l 
[1  + b^(27T  kL)^] 


where 


n = 1/2  (-m-1)  and  b = 

The  use  of  these  equations  first  implies  knowledge  of  m and  L.  These 
parameters  displayed  considerable  variation  between  individual  turbulence 
runs  in  the  MEDCAT  data.  However,  for  the  purposes  of  presenting  the 
rms  data,  it  was  assumed  that  L = 3000  feet  and  m = -1.83  for  the  lateral 
component  and  m = -1.33  for  the  vertical  component. 

2 

Histograms  of  the  estimates  of  o are  presented  in  Figure  18.  In 
interpreting  these  higtograms  it  is  important  to  realize  that  the 
editing  criteria  for  the  turbulence  runs  of  this  figure  included  the 
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requirements  that  > 1.5.  Further,  in  collection  of  the  data,  the  pilot 
sought  the  atmospheric  regions  predicted  as  being  most  likely  to  contain 
clear  air  turbulence.  Thus,  this  histogram  cannot  be  interpreted  as  a 
random  sample  of  rms  values  since  the  data  base  is  not  a random  sample 
and  the  truncation  point  of  the  o distribution  was  variable. 


SECTION  V 
PROJECT  HICAT 


The  HICAT  program  results  are  based  upon  flight  measurements  of 
aircraft  eg  acceleration  response  and  the  true  gust  velocity  components 
together  with  measurements  of  the  related  aerodynamic  and  meteorological 
variables.  HICAT  tests  involve  deliberate  flights  into  regions  of 
turbulence  or  suspected  turbulence.  A major  objective  was  to  determine 
gust  velocity  spectra  for  turbulence  ranging  from  100  feet  to  50,000  feet 
in  wave  length.  These  flights  were  designed  to  sample  the  variety  of 
atmospheric  phenomena  associated  with  turbulence  in  the  40,000  to  70,000  ft 
altitude  region.  This  program  provides  the  basic  data  for  the  development 
of  a statistical  model  of  world-wide  high-altitude  clear-air  turbulence. 

The  following  tabulation  from  Reference  9 summarizes  all  the  HICAT  flight 
experience  including  that  from  the  initial  and  redirected  HICAT  programs; 


HICAT 

Program 

Test  Period 

Tests 

Tests 

Encountering 

HICAT 

FI ight 
Hours 

HICAT 

Hours 

Initial 

Feb  64  - 

Oct 

64 

33 

17 

94.5 

7.4 

Redirected 

Oct  65  - 

Feb 

67 

146 

84 

649.5 

29.2 

Extended 

Mar  67  - 

Feb 

68 

106 

72 

477.6 

18.3 

All 

Feb  64  - 

Feb 

68 

285 

173 

1221.6 

54.9 

An  extensive  analysis  of  the  meteorological  and  geophysical  conditions 
with  high  altitude  clear  air  turbulence  was  conducted  and  reported  in 
Reference  10.  The  distributions  of  the  number  of  observations  by  season 
and  latitude  are  shown  in  Figure  19.  The  ratio  of  flight  miles  in 
turbulence  to  the  total  flight  miles  was  analyzed  as  a function  of  topography. 
Four  categories  were  defined  as  follows: 

Water 

Flatland  (local  relief  differences  • 3000  ft.) 

Low  mountains  (local  relief  differences  3000  to  7000  ft.) 

High  mountains  (local  relief  differences  ■ 7000  ft.) 
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Figure  19.  Distribution  by  Season  and  Altitude  of  Number  of  Observations 
of  High  Altitude  Clear  Air  Turbulence 
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Figure  20  shows  the  distribution  of  the  number  of  observations  by  topography 
and  by  season.  HICAT  data  indicate  that  the  percentage  of  turbulence  is 
related  to  the  topography  immediately  below  the  aircraft.  Figure  21 
presents  the  distribution  of  total  flight  miles  and  turbulent  flight 
miles  by  terrain  for  ferry  flights  and  all  flights. 

An  important  objective  of  the  HICAT  program  was  to  derive  gust 
velocity  spectra  for  turbulence  waves  ranging  from  100  feet  to  as  much 
as  50,000  feet  in  length.  To  achieve  uniform  statistical  reliability  it 
is  required  that  sample  lengths  increase  in  proportion  to  the  largest 
wavelength  of  interest  in  the  spectrum.  The  square  root  of  the  integral 
of  the  spectral  curve  from  the  short  wavelength  limit  to  the  longest 
wavelength  of  interest  for  the  design  and  operation  of  aircraft  is  an 
important  quantity,  but  one  that  is  rarely  obtainable  from  direct  measure- 
ments. A method  was  presented,  in  References  9 and  11  for  comparing 
spectra  with  different  limits  using  the  root  mean  square  (rms)  values  for 
the  following  wave  length  limits:  1000,  2000,  4000,  10,000,  20,000,  and 

40,000  feet.  The  short  wavelength  limit  was  established  at  140  feet. 

The  envelopes  of  minimum  wavelengths  and  maximum  wavelengths  for  turbulence 
investigations  are  presented  in  Figures  22  and  23  from  Reference  12. 

The  ratios  of  the  rms  at  A = 2000  ft,  A = 4000  ft,  A = 10,000  ft, 
and  A = 20,000  ft  to  rms  at  A = 1000  ft  were  used  as  quantitative 
measures  of  the  shape  of  the  spectra.  Table  6a  presents  the  average 
values  of  these  ratios  for  two  categories  of  topography  and  for  three 
values  of  the  maximum  wavelength  of  the  spectra.  Table  6b  also  includes 
the  calculated  ratios  for  the  "mild  knee"  equation  with  the  slope  m of  the 
short  wavelength  equal  to  -5/3. 

The  "mild  knee"  equation  is 

i'iQ)  = 

1 + (fiL)'"' 

when  T’  { ) is  the  power  spectral  function,  J is  the  reduced  frequency 
(2n/  ■)  and  L is  the  scale  of  turbulence.  This  equation  appeared  to  best 
fit  the  HICAT  spectra  data  for  the  three  components  of  the  gust  velocity 
presented  in  Table  6a. 
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Figure  20.  Distribution  of  the  Number  of  Observations  of  Turbulence  by 
Topography  and  by  Season 
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IL 


FERRY  FLIGHTS 


ALL  FLIGHTS 

Figure  21.  Distribution  of  Total  Flight  Miles  and  Turbulence  Flight 
Miles  by  Terrain  for  Ferry  Flights  and  All  Flights 
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TABLE  6a 


AVERAGE  VALUES  OF  THE  RATIO  OF  U^,  Up,  AND  Up 
AT  RMS  (2000),  RMS  (4000),  RMS  (10,000),  AND  RMS  (20,000)  TO  RMS  (1000) 


X max 
4000  ft 


X max 
10,000  ft 


Samples 


RMS  (2000) 
RMS  ClOOOi 

RMS  (4000) 
RMSTTOOO) 


Sampl es 

RMS  (2000) 
RMS  (1000) 


RMS  ±4000] 
RMS  (1000) 


RMS  (10,000J 
RMS  (1000) 


Samples 


Water,  Flatland 

^'v 

Uf 

^L 

(13) 

(31) 

(44) 

X max 
20,000  ft 


RMS  1 

BBSS 

1 

RMS  1 

1 

RMS  1 
RMS  1 

4000 ; 
: 1 000 : 

) 

i 

RMS  1 
RMS  ( 

;io,ooo) 
: 1 oooT 

RMS  (20,000) 
RMS  OOOO') 


± ± ‘^L 


1.26  I 1.25  I 1.29  1.32  1.30 

1.65  1.69  1.77 


(16)  (16)  (19)  (20) 

1.25  1.29  1.32  1.28 

1.54  1.63  1.67  1.75 


1 .65  1.95  1 .91  2.06  2.31  2.45 


(4)  (4) 


(5)  (5) 


1 .21  I 1 .29  1.24  1 .31  1.34  1.29 

1.36  I 1.53  1.50  1.72  1.75  1.76 


1.54  1.80  1.71  2.25  2.51  2.50 

.67  2.01  1.86  2.66  3.19  3.11 


A total  of  732  samples  for  which  the  fliqht  path  was  level  and  th( 
turbulence  was  indicated  to  be  very  light  were  used  to  determine  the 
distribution  lengths  of  the  turbulent  regions.  The  percentages  of  ex- 
ceedances were  computed  for  the  total  sample  for  cases  of  moderate  or 
greater  turbulence  (144  cases)  and  severe  turl)ulence  (16  cases). 
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THEORETICAL  VALUES  OF  THE  RATIOS  OF  RMS  (2000), 

RMS  (4000),  RMS  (10,000)  AND  RMS  (20,000)  TO  RMS  (1000) 
FOR  VARIOUS  L'S  OF  THE  MILO  KNEE  EQUATION  WITH  m - -5/3 


RMS  (10,000) 
RMS  (1000) 


RMS  (20,000J 
RMS  (1000) 


500  1000  2000  4000 


1.23  1.26  1.30  1.32 


1.60  1.67 


.96  2.17 


1.68  1.92  2.22  2.55  2.58 


The  mean  lengths  of  the  turbulent  regions  by  altitude  interval  and  season 
are  given  in  Tables  7 and  8.  The  data  from  these  tables  indicate  that 
the  lengths  of  the  turbulent  regions  decrease  with  increasing  altitude 
and  that  the  lengths  of  the  turbulent  regions  are  greater  in  winter  than 
in  the  other  three  seasons. 


TABLE  7 


AFFDL-TR-77-4 


TABLE  8 


MEAN  LENGTHS  OF  TURBULENT  REGIONS  AS  A 

FUNCTION  OF  SEASON 

Season 

No.  of  Cases 

Mean  Length 

Winter 

219 

29.5  nm 

Spring 

196 

22.5  nm 

Summer 

81 

22.5  nm 

Autumn 

236 

21  . 5 nm  j 

In  Reference  13,  the  numerical  values  for  the  standard  deviation  of 
the  rms  gust  velocities  (b-|),  the  scale  lengths  and  slopes  of  the  power 
spectral  density  curves,  and  the  ratio  of  turbulent  flight  miles  to  total 
flight  miles  were  determined.  These  numerical  values  are  given  in  Table  9. 
The  values  of  b^  are  shown  to  increase  significantly  with  roughness  of 
terrain  for  all  three  components.  Numerical  values  are  0.9  to  3.4  feet 
per  second  for  the  vertical  component,  1.5  to  4.6  feet  per  second  for  the 
lateral  component,  and  1.2  and  4.6  feet  per  second  for  the  longitudinal 
component.  Table  10  lists  the  ratios  of  the  turbulent  flight  miles  to 
total  flight  miles  that  are  recommended  for  use  in  the  design  and  operation 
of  advanced  aircraft.  The  ratios  obtained  from  an  analysis  of  the  HICAT 
data  were  judged  to  be  high  mainly  because  of  the  search  procedures  used. 
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TABLE  10 

RECOMMENDED  RATIOS  OF  TURBULENT  FLIGHT  MILES  TO  TOTAL  FLIGHT  MILES 


Category 

Ratio  Turbulent 

65,000-70,000  Ft 

0.010 

60,000-65,000 

0.010 

55,000-60,000 

0.020 

50,000-55,000 

0.030 

45,000-50,000 

0.030 

45,000-70,000 

0.020 

Water 

0.018 

Flatland 

0.018 

Low  Mountains 

0.022 

High  Mountains 

0.030 

Winter 

0.024 

Spring 

0.024 

Summer 

0.014 

Autumn 

0.018 

I 


50 


AFFDL-TR-77-4 


SECTION  VI 
CONCLUSIONS 

The  turbulence  data  obtained  from  the  ALLCAT  Program  have  been 

integrated  with  data  from  all  other  available  sources  to  develop  improved 

gust  design  criteria.  Gust  design  procedures,  based  on  power  spectral 

methods  of  generalized  harmonic  analysis,  reguire  modeling  of  the 

atmospheric  turbulence  environment  in  terms  of:  a spectrum  shape,  a 

root-mean-square  gust  intensity  (o  ),  a scale  of  turbulence  (L),  and  the 

w 

parameter  (P)  depicting  the  proportion  of  flight  time  spent  in  turbulence 
(see  References  14  and  15).  Updated  P and  o values  are  presented  in 
Figures  24  and  25.  The  aircraft  response  to  turbulence  is  treated  in 
the  frequency  plane  as  shown  by  the  equation 

= i H^(a)) 

which  relates  the  output  spectrum  of  the  response  variable  x to  the  input 
spectrum  of  gust  velocities  through  the  frequency  response  function 
defined  as  the  response  x that  results  from  a sinusoidal  gust  encounter 
of  unit  amplitude.  This  response  equation  leads  to  the  two  basic  struc- 
tural response  quantities  A and  N^,  where  A relates  the  root-mean-square 
output  value  o to  the  root-mean-sguare  gust  input  o,,  value  throuah  the 
relation 


and  N^  defines  the  number  of  zero  crossings  of  x with  positive  slope  per 
second. 

Estimates  for  the  turbulence  scale  L have  ranoed  from  200  ft  to 
over  5,000  ft.  The  effective  scale  of  turbulence  for  all  altitudes  is 
recommended  as  L = 750  ft  (see  Reference  14).  The  recommended  spectrum 
for  the  vertical  turbulence  is  the  von  Karman  sp'ctrum. 

Load  exceedance  curves  have  been  established  for  design  based  on 
mission  considerations.  These  updated  curves  are  presented  in  Figure  26. 
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Figure  26.  (Continued) 
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